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Exploring the fundamental quantum behaviour of optomechanical resonators is of great interest recently but 
requires the realization of the strong coupling regime. We study the optical photon-phonon coupling of the 
so-called membrane in the middle (MITM) optomechanical system. Using coupled-mode theory we find that 
the optomechanical coupling is proportional to the electric susceptibility of the membrane. By considering 
the doping atoms or spins into the membrane and driving these appropriately we induce a tunable ultra-large 
refractive index without absorption which enhances the optomechanical coupling. Using this we predict an 
ultra-strong single-optical photon strong coupling with large quantum cooperativity for Er 3+ dopants at low 
temperature, while Cr 3+ in a Ruby membrane may display ultra-large quantum cooperativity at room tempera¬ 
ture. Our scheme also can tune the strength of the coupling over a wide range and can also control whether the 
optomechanical force is attractive or repulsive. Our work opens a door for fundamental physics and applications 
relying on the realization of the strong coupling regime in quantum optomechanical systems. 

PACS numbers: 85.60.Gz, 42.50.Dv, 42.65.-k 


Optomechanical systems have made tremendous progress 
over the past decade. It is a promising testbed for fundamental 
quantum mechanics from manipulating the motional quantum 
ground state of mescoscopic objects JT] [2], to non-Gaussian 
quantum state synthesis 0 , and also a platform for important 
applications such as squeezed light a, ultrasensitive mea¬ 
surement a, diabolical points 0, and microwave(mw)-to- 
optical quantum interface and quantum networks 00. 

The above schemes crucially rely on amplifying the rather 
weak single-photon coupling g om . Observation of a non- 
Gaussian state of the mechanical resonator even requires 
a single-photon coupling stronger than the linewidth k of 
the cavity, i.e. the so-called single-photon strong coupling 
(SPSC) regime J3}. SPSC has only been demonstrated so 
far in optomechanical arrangements utilizing ultracold atomic 
gases 0. SPSC continues to remain a very challenging 
regime for mescoscopic mechanical resonators. The state-of- 
the-art current experiments in mesocopic optomechanics typ¬ 
ically can only achieve the small coupling-decay ratio (CDR) 
of 1(T 3 Q0). Even after amplification by a large optical coher¬ 
ent state | or), the record experimental CDR for optical setups 
~ i.5 mi- On the other hand, there is much interest in the 
quantum behavior of optomechanical systems which have a 
quantum cooperativity (QC) in excess of unity BT0HT21 . In this 
case the light and mechanics can coherently exchange faster 
than the thermal decoherence of mechanical motion. Only 
very recently, two theoretical proposals making use of ultra¬ 
sensitive superconducting circuits ltl3l|T4l approach the SPSC 
regime for microwave (mw) photons, g om > k, but require 
cryogenic temperatures. An array of identical mechanical res¬ 
onators supporting a collective oscillation is also theoretically 
studied for strong coupling to an optical mode BE), while a 
very recent theory proposal aims to achieve SPSC via an op¬ 
tically trapped nanodiamond ED- However, in typical setups 
the coupling between a single mechanical mode and optical 
mode is weak. 

In this letter, we propose a dynamically tunable optome- 




FIG. 1. (Color online), (a) Schematic of the optomechanical res¬ 
onator. A membrane is inserted into an optical Fabry-Perot cavity. 
An ensemble of A-type atomic systems are doped in the membrane, 
(b) The A-atomic configuration is optically but incoherently pumped 
from its ground state |1) to the excited state |3). A coherent mi¬ 
crowave field detuned by A,, is applied to create coherence between 
the doublet ground states |1) and |2). The probe (cavity mode a) field 
detuned by A,, drives the transition |2) «-> |3). 


chanical system, see Fig. [I] We find that the optomechanical 
coupling is proportional to the electric susceptibility x of the 
membrane. Based on this we propose a theoretical scheme 
to enhance the electric susceptibility of the mechanical res¬ 
onator by several orders in magnitude. We greatly enhance 
the optomechanical coupling without introducing extra loss to 
the cavity, to achieve the SPSC regime. Ultrarefractive media 
have been studied in various systems GM2, and has been 
demonstrated in experiment Ii20ll . but it has not been exploited 
for enhancing optomechanical coupling so far. 

One implementation of our scheme is depicted in Fig. [lja). 
A membrane with thickness / is inserted at equilibrium posi¬ 
tion zo in a Fabry-Perot(FP) cavity with length L and finesse 
F , a setup well studied in experiments lf2Tl . The membrane 
has mass M and oscillates with frequency Q„,, with zero-point 
motion z zp = yJfi/2MQ. m , and we represent the quantised 
membrane’s displacement from zo as z m = z zp (fc l + b ), where 
bHb) are the creation(annihilation) operators of the mechani¬ 
cal motion. The membrane is made from a host medium with 
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a relative permittivity s h. An ensemble of A-type quantum 
systems is doped in the membrane to create extra susceptibil¬ 
ity x, and whose atomic level system is schematically shown 
in Fig. |T]b). This configuration can induce a giant electric sus¬ 
ceptibility to substantially modify the phase of a probe field 
passing through it. 

Mode coupling Using coupled-mode theory |22j 123], we 
present a derivation of the optomechanical coupling strength 
involving many more parameters and from this we will ob¬ 
serve the dependence on the electric susceptibility of the 
membrane. We first derive the coupling between different cav¬ 
ity modes without the mechanical motion of the membrane. 
We consider the membrane in Fig. |TJa) as a slab of medium 
with a macroscopic complex polarization P - sq(x+si, -1 )E p , 
where £q is the vacuum permittivity, and E p is the probe field 
input into the membrane. The overall relative permittivity is 
X + £/,. We set £ - £' + £”, where £' = Re [<(-'] and £” = Im[£]. 
We consider now an empty optical FP cavity. The total elec¬ 
tric field trapped in the cavity is E p , and we seperate this into 
the positive and negative frequency components, E = E' +E , 
£ ,± = e T; " f £(r), where E( r) is the E-field spatial distribution 
at position r in the cavity. This E-field obeys Maxwell’s wave 
equation V 2 E - 4 ^ 8 = 0, where e,.(r) is the overall rel¬ 

ative permittivity of the optical cavity at position r, c is the 
speed of light in vacuum and e r ( r) = 1, in the absence of the 
membrane. The E-field in the cavity can be decomposed as a 
superposition of eigenmodes as 

E + = di(t) E + ( r, f) = £ 3 kit) e- j0J " E,{ r) , (1) 

/ / 


where E/( r) is the E-field distribution for the Vth eigen¬ 
mode, a>i is the resonant frequency and <?/ is the polariza¬ 
tion of this mode and ) is the annihilation and creation 
operator of the /’th mode. For each eigenmode, we have 
c 2 V 2 EJ-d 2 E{)E^/8t 2 = 0. Substituting Eq. |lj into Maxwell’s 

i ^ c , 1 d 2 a/£ r EUr,t) „ 

wave equation, we have 2_,/ a / V EJ(r,t) - - -^4 - = 0. 

It is also reasonable to assume that s r is time-independent in 
the absence of a large time-dependent driving of the mem¬ 
brane. Normally, we have |-2yo>/^| » \ leading to 


Thus, Maxwell’s wave 
equation, as a decomposition of the eigenmodes, has the form 


\ ds r ^{r,i) dai „ d 2 (s h - 1 + x)E J [(r,t) 

A | 2 J, aT + »'-*3- 


We focus on the modification of the Ath mode by the 
mechanical motion but neglect the intermode cou¬ 
plings in this paper. Using the orthogonality condi¬ 
tions f e k ■ ei El(r)Ei(r)dr = 6 k i, the coupled mode 
equation in the absence of mechanical motion reads as 
= jAa> k d k - AK k a k , where A 64 and A K k corresponds 
to the frequency shift and linewidth change by insert¬ 
ing the membrane, respectively. One finds 2U k Ato k — 
o) k e Q j v ^x'\E k (r)\ 2 dY + oj k s 0 / (e^ - l)\E k (r)\ 2 dr and 


2U k Ax k = u> k so (s'! + x")\E k (r)\ 2 dr, with the volume 

^Fobj n 

of the membrane V 0 bj and U k = so C, \E k (r)\ 2 dr. These 
formula are in agreement with those derived for toroidal 
cavities 1231 . Defining the E-field E k (r) = E k $e k (r), 
where the amplitude E k _o = y/tiu) k /soV k , with the mode 
volume V*, and e k (r), the spatial distribution of the mode 
(with f \e k (r)\ 2 dr/V k = 1), we take the normalized 

mode distribution to be e k (r) = f k (x, y, z) sm(k k z), where 
f k (x,y,z ) is the transverse profile of the k’th mode, vary¬ 
ing slowly along the z-axis, and k k is the z-component of 
the propagation constant. Here we assume that the field 
propagates along the z-axis such that k k = 2n/A k with the 
wavelength A k = 2nc/u) k . In the FP cavity, sin(k k L) = 0. 
In the case of l <s A k , f k (x,y,z ) ~ f k (x,y,zo) in the mem¬ 
brane. We have U k = ^so\E k o\ 2 A k with the effective mode 
cross section area A k = 2 V k /L. For some position z 0 , 
A k = f f \f k (x, y, zo)\ 2 dxdy. 

Now we consider the mechanical motion of the mem¬ 
brane. The displacement d m of the membrane around 
its equilibrium position zo is dependent on the mechanical 
mode. The displacement can be described as d m (x,y,z ) = 
Z m fm(x,y), where f„(x,y) is the transverse mode distribu¬ 
tion of the mechanical vibration and z m is the displacement 
amplitude. Correspondingly, the electric field of the kth 
cavity mode at the displaced membrane location becomes 
e k (r) = f k (x,y,z 0 + d m ) sin k k (z 0 + d m ) ~ f k (x,y,z 0 )[s'mk k z 0 + 
k k cos k k zofm(x,y)z m ], expanding to first order in d m . 

Applying z m = z /p (b' +b), the optomechanical Hamiltonian 
for a specific cavity mode a takes the form, when kl « 1, 

H m = [u) c + Ao> + i(A k + /c)] a 1 a 

+ [gom,p + gom.h] a f a(b' + b) , 

with k is the total decay rate of cavity without the membrane. 


Aco « (j) c l(s' h - 1 +x') sin 2 (fco)/T , (4a) 

Ak * Q c kI (s'l + x") sin 2 (fe 0 )/L , (4b) 

gom.p = X^c Zzp kl sin(2fe 0 )S/L , (4c) 

<?om,h — (s /? - 1) 0) c z z p kl sin(2kzo)S/L , (4d) 


where S — f f \f k (x,y,zo)\ 2 fm(x,y)dxdy/A k can be almost 
unity if f„(x,y) = 1 within the optical mode waist, and 
Q c = o>Jk is the quality factor of the cavity. According to Eq. 
0 . the coupling is not only determined by ojJL but also pro¬ 
portional to the electric susceptibility of the material. There¬ 
fore, a giant susceptibility x' without extra additional absorp¬ 
tion, i.e. x" - 0 ’ can greatly enhance the optomechanical 
coupling strength. 

We choose to separate the couplings into a contributions 
due to the host material, g om ,h, and the induced polarization 
g om p due to the doped atomic systems. The total optome¬ 
chanical coupling is g om = g om ,p + gom.h- We now explore 
how our scheme improves the of-interest coupling-decay ra¬ 
tio (CDR). We assume that the intrinsic loss rate of the cavity 
is Kj. With critical external coupling we have K ex = Kj and 
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k = /c ex + Kt, and in the absence of the atomic polarization the 
CDR is gom,h/ K - Considering now the induced polarization of 
the atoms the intrinsic decay becomes k' = a - ,- + A k and the 
total decay changes to k' = 2 k, + A k under the critical cou¬ 
pling condition. After including the change in decay rates, the 
improved CDR becomes 

§om _ gom.h 

k' K'( e h - 1 ) K 

We are interested when, \x\ » (s/,-1), yielding G ~ xl( e h~ !)• 
Giant susceptibility We now present a method to induce 
a giant polarization without extra absorption (and potentially 
gain), using a A-type atomic level system, see Fig. |TJb). The 
advantages are twofold: (i) the giant susceptibility^ can en¬ 
hance the optomechanical coupling; (ii) any gain can reduce 
the decay rate of the cavity. For this purpose, a small popu¬ 
lation in the excited state and a quantum coherence between 
the ground states are essential Q3. To proceed we assume 
the system is initially populated in |1) but is incoherently 
pumped to the excited state | 3 ) with rate r. A mw field Q /( 
is applied to create coherence between the two ground states. 
The cavity mode a is modeled as a probe field, Q ; , = QE p /Ti 
with E p = Eo(a), where Q is the dipole moment of tran¬ 
sition |2) <-» | 3 ). For simplicity, we assume the atoms are 
identically polarized so that the complex macroscopic polar¬ 
ization is P = NQf>r\ = S\)X V E V , where N is the number 
density of atoms, p is the density matrix and p,j = (i\p\j), 
and Xp denotes the susceptibility due to the atoms as seen 
by the probe field. Inducing a large Xp requires /V/l 3 » 1. 
In such dense dielectric medium the macroscopic probe field 
E p that couples to atoms must be replaced by the local mi¬ 
croscopic electric field El, which is related to the macro¬ 
scopic volume polarization P by the Lorentz-Lorenz rela¬ 
tion fl8l [24], namely E/ = E ; , + P/ 3 eo- This near dipole- 
dipole(NDD) effect leads to an effective susceptibility for E p , 
Andd = CD- The time evolution for the macroscopic 

density matrix is determined by p - -i(Hp-pH)+X?p, where 
H = -A p a-22-(Ap+Ap)o-ii+£2p(cr32+o-23)+Qp(cr2i+o-i2) with 
operators cr,- ; = |i)(y|, and SA’p = y/2{2ApA 4 - A^ Ap - pA f A} 
with y = {71,72, r} corresponding to A = {o"i 3 ,cr 2 3 ,cr 31 }. 
71(72) is the decay rate from | 3 ) to |1)(|2)). The decay rate 
can be estimated as 7 = 4w 3 <3 2 y/si,/6n:eobc 3 . For simplicity, 
we assume 71 = 72 = 7 and neglect the decoherence of |2). 
For a first-order analysis in the probe field E p liTTl |19l , we 
obtain the macroscopic polarization in the steady state 


-w) 

[2/(r-2/A„) 

[r + 7 - 2i(A p + A^)] 

1 - 

• 8 / 0 2 } 

(r - 2iA p ) | 

[(7 - 2iA p ) 

[r + y- 2i(A p + A^)J 

| + 4^} 


where so = NQ 2 /sotiy = 3/VT 3 /8;r 2 y/sh ~ NA 3 /26 yjsh. 

To obtain a large increase in the optomechanical coupling 
we need to accurately control the doping of the atoms so that 
[^ p /3 - 1 | is small but much larger than \x p V which then yields 
Tndd ~ A p /(1 -T p /3) »A P - 

Replacing^ with ^ndd in Eqs. (|4]i and ([5]), we find that the 
coupling gom.p Andd, and thus a giant value of [^ndd[ is ob- 




FIG. 2. (Color online). Real (blue solid lines) and imaginary (red 
dashed lines) parts of the electric susceptibility ^ NDD as a function of 
the detuning A p (a) Cl p /y = 1 and (b) Cl p /y = 0.1. s 0 = 3/1.66, r = 
0.1 7 , = 0.4y for^NDD- x" < 0 implies gain. 



Rp/r Rp/y 


FIG. 3. (Color online). Enhancement of optomechanical coupling- 
decay ratio (CDR): (a) |g om //c'| as a function of the pump rate r and 
the mw field Here so = 3/1.66, A p = 0.3y, A p = 0.4y. (b) 
Enlarged central area of (a). Contours (dashed) are guide to the eye. 
In (b) (or (a)) dashed gray circle, overlapping with white spot, shows 
\goml K '\ > 1 or 0.1. Region above(below) the black line is gain(loss) 
and the yellow contour shows significant loss W^x" > 0.1 and gain, 
green contour 10 ~ 3 ;y" < - 0 . 1 , regions. 


tained. The argument of the complex ^ndd 6 C, can be used 
to control the sign of the coupling, g om . In Fig. [2] we plot 
Andd- for both large and small mw driving, Q /( /y = { 1 , 0 . 1 }. 
In the case of large mw driving, and A p < 0.303y, we find 
significant enhancement [TnddI > w f | ^ c th e medium is 
transparent or displays gain for the probe field (Tndd — 0 )- 
This induced gain can considerably reduce the intrinsic loss 
of the cavity and subsequently leads to k' <s k lf25l |26| . As 
a result, the CDR can be improved by more than three or¬ 
ders if (e/, - 1) < k/k' , see Eq. d5 1 . If a weaker mw drive is 
applied, e.g. El p /y = 0.1, (see Fig. 2lb)), the atomic suscepti¬ 
bility drops precipitously and we switch off the optomechan¬ 
ical coupling. Here sq = 3/1.66 requires a number density 
N = 10 14 - 10 15 cm~ 3 , for A = 0.5 - 1 pm, corresponding to 
AT A 3 -100. 

We can now estimate the possible CDR that can be achieved 
by our scheme, g 0 mli <'. Typically, in normal media g 0 m,h/^ ~ 
10 3 and Eh — 1 < 10 , but due to the induced gain in our 
scheme, k/kAeu - 1) > 10 is possible. Thus it is reasonable 
to assume = ICC 3 considering only the gain. When we 

fix the detuning, A p = 0.3 7 , A p = 0.4y, and the concentration 
of atoms, so = 3/1.66, the SPSC regime, gem/*' > 1, can be 
achieved over a large range of pump and mw drivings while 
the medium displays gain, see Fig. [3] Normally, the density 
of atoms is fixed once the device is made and is therefore hard 
to control accurately. Figure |4j a) shows that SPSC can be re¬ 
alized over a range of number densities 1.583 < so < 1-625, 
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FIG. 4. (Color online). Enhancement of the CDR over a range of 
concentration for the atomic dopant in the resonator. Modulus (a) 
and (b) argument [scaled by tt], of the ratio g om /K' for different mw 
driving rQ p and number density N. We choose r = 0.098y, A p = 
0.3y, A p = 0.4y. Black contour shows the strong coupling region 
gom/K' > L while green contour for a gain medium, W^x" < -0.1 
(The gain inside the contour is larger). 


by tuning the mw driving (overlap region between green and 
black contours). Interestingly, the phase of the coupling g om , 
switches from 0 <-> n, indicating a change in sign of the op¬ 
tomechanical coupling when the mw driving is tuned across 
the white-blue boundary in Fig. |4jb), due a change in sign 
in _y' NDD . Note that the direction of the radiation pressure, 
F rp = -hg mn (a'a)/z zp is dependent on the sign of the cou¬ 
pling |[27l . and we find therefore, that the radiation pressure 
can be dynamically tuned from being repulsive to attractive. 

Now we present possible implementations using existing 
experimental systems. Atomic A-type systems have been 
demonstrated in various solid-state systems from quantum 
dot m . Nitrogen-vacancy (NV) |29||30l or Silicon-vacancy 
(SiV) ||3T1 . in diamond to atoms l(32l . and rare earth ions l33l . 
For our quantum configuration, we consider Er 3+ ll25l {331 
at low temperature (T = 10 K), implanted in Silicon Ni¬ 
tride (Si 3 N 4 ), membrane with Sh - 4 |34j. By choosing 
so = 3/1.66 corresponding to N « 2.57 x 10 13 cm -3 , and 
A p = 0.30285y, at wavelength A = 1550 nm, we can induce a 
susceptibility of Undd = 1181.45 - 0.70/, at sin 2 (feo) = 0.5, 
when r jy = 0.1, and D /( /y = 1. For the cavity, we take 
L = lOOT, F — 2 x 10 5 , yielding Q c - 2 x 10 7 (/e/27r = 
9.7 MHz),, and k/k' — 30 l35l For a Si 3 N 4 membrane 
with thickness / = 100 nm, diameter I) - 10 pm, mass 
density p m — 2.7 g-crrT 3 (m = 21 pg), a tensile stress of 
T s ~ 0.9 GP and Q m = 4x 10 6 , ED El E3, yields a mo¬ 
tional oscillation frequency Cl m /2n = 40.8 MHz, zero point 
motion z zp = 3.1 fm, and B = 0.92. From this we obtain a 
very large CDR, g Qm /x' = 5.3, and more importantly, an ex- 
tremly large quantum cooperativity Cq = gom/^Tm = 174.7 
at the single-photon level, where y/, = y,„ x n± is the me¬ 
chanical decoherence rate at 10 K = 5.1 x 10 4 ) MH3. 
If we take Cr 3+ (A = 694 nm) in Ruby (/ = 50 nm,p,„ = 
3.98 g ■ cm -3 , T s = 0.3 GPa) at T - 300 K f33l but assume 
lower finesse of F — 2 x 10 4 yielding Q c = 2 x 10 6 as an 
estimation of absorption in Ruby ll37ll . we can still achieve 
gom/K 1 = 22.4 and Cq ~ 230.8 at room temperature. As men¬ 
tioned above, Cq > 1 allows a number of coherent optical 
control techniques to manipulate the quantum state of the me¬ 


chanics. 

In conclusion, based on our coupled mode analysis a giant 
enhancement of the optomechanical coupling can be achieved 
by inducing an ultrarefractive index in the membrane in the 
middle based optomechanical system. This giant enhance¬ 
ment permits low temperature explorations of the fascinat¬ 
ing regimes of SPSC and QC in excess of unity and for cer¬ 
tain materials, e.g. Ruby, a room temperature implementation 
may be possible. By tuning the pump laser or mw driving we 
can rapidly switch on and off the optomechanical coupling 
or change between repulsive and attractive optomechanical 
forces. 


* keyu.xia@mq.edu.au 
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